Particle emission in hydrodynamics: a problem needing a solution 



F. Grassi 

Instituto de Fisica, Universidade de Sao Paulo, C. P. 66318, 05315-970 Sao Paulo-SP, Brazil 

A survey of various mechanisms for particle emission in hydrodynamics is presented. First, in the 
case of sudden freeze out, the problem of negative contributions in the Cooper-Frye formula and 
ways out are presented. Then the separate chemical and thermal freeze out scenario is described 
and the necessity of its inclusion in a hydrodynamical code is discussed. Finally, we show how 
to formulate continuous particle emission in hydrodynamics and discuss extensively its consistency 
with data. We point out in various cases that the interpretation of data is quite influenced by the 
choice of the particle emission mechanism. 
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1. INTRODUCTION 



Historically, the hydrodynamical model was suggested 
in 1953 by Landau pj] as a way to improve Fermi statis- 
tical model 2]. For decades, hydrodynamics was used to 
describe collisions involving elementary particles and nu- 
clei. But it really got wider acceptation with the advent 
of relativistic (truly) heavy ion collisions, due to the large 
number of particles created and its success in reproducing 
data. Brazil has a good tradition with hydrodynamics. 
Many aspects of it have been treated by various persons. 
For illustration, the following papers can be quoted. Ini- 
tial conditions were studied in 0|. Solutions of the 
hydrodynamical equations using symmetries 5, 6] or nu- 
merical 0, m were investigated. The equation of dense 
matter was derived in |g|.ll(J. Comparison with data was 
performed in 0, 0, ll3L Il4 Il k 16, 1 71. The emission 
mechanism was considered iiTQ, Ell El Eaf. 

In this paper, I concentrate on the problem of particle 
emission in hydrodynamics. In the Fermi description, en- 
ergy is stored in a small volume, particles are produced 
according to the laws of statistical equilibrium at the 
instant of equilibrium and they immediately stop inter- 
acting, i.e. they freeze out. Landau took up these ideas: 
energy is stored in a small volume, particles are produced 
according to the laws of statistical equilibrium at the in- 
stant of equilibrium, expansion occurs (modifying parti- 
cle numbers in agreement with the laws of conservation) 
and stops when the mean free path becomes of order the 
linear dimension of the system, which led to a decoupling 
temperature of order the pion mass for a certain energy 
and slowly decreasing with increasing energy. In today's 
hydrodynamical description, two Lorentz contracted nu- 
clei collide. Complex processes take place in the initial 
stage leading to a state of thermalized hot dense matter 
at some proper time tq. This matter evolves according 
to the laws of hydrodynamics. As the expansion pro- 
ceeds, the fluid becomes cooler and more diluted until 
interactions stop and particles free-stream towards the 
detectors. In the following, I review various possible de- 
scriptions for this last stage of the hydrodynamical de- 
scription. The usual mechanism for particle emission in 
hydrodynamics is sudden freeze out so I will use it as a 



point of comparison. I will start in section 2, reminding 
what it is, some of its problems and ways outs. There 
is another particle emission scenario which is a small ex- 
tension of this idea of sudden freeze out: the separate 
chemical and thermal freeze out scenario. It has become 
used a lot e.g. to analyse data. So I will discuss in section 
3 what it is, its alternatives and how to incorporate it in 
hydrodynamics. Continuous emission is a mechanism for 
particle emission that we proposed some years ago. As 
the very name suggests, it is not "sudden" like the usual 
freeze out mechanism. I will explain what it is precisely 
in section 4 and how it describes data compared to freeze 
out. Finally I will conclude in section 5. 



2. SUDDEN FREEZE OUT 
2.1 The traditional approach and its problems 

Traditionally in hydrodynamics, the following simple 
picture is used. Matter expands until a certain dilu- 
tion criterion is satisfied. Often the criterion used is 
that a certain temperature has been reached, typically 
around 140 MeV in the spirit of Landau's case. In 
some more modern version such as |2^ . a certain freeze 
out density must be reached. There also exist attempts 
|l9l I25L l2a. I27I l2a | to incorporate more physical informa- 
tions about the freeze out, for example type i particles 
stop interacting when their average time between interac- 
tions rj^^jj becomes greater than the fluid expansion time 
and average time to reach the border. When the freeze 
out criterion is reached, it is assumed that all particles 
stop interacting suddenly (this is called "freeze out" ) and 
fly freely towards the detectors. As a consequence, ob- 
servables only reflect the conditions (temperature, chem- 
ical potential, fluid velocity) met by matter late in its 
evolution. 

In the sudden freeze out model, to actually compute 
particle spectra and get predictions for the observables, 
the Cooper-Frye formula ^ |Q may be used. 



E(fN 



(1) 



2 



dcr^ is the normal vector to this surface, the particle 
momentum and / its distribution function. Usually one 
assumes a Bose-Einstein or Fermi-Dirac distribution for 
this /. 

This sudden freeze out approach is often used however 
it is known to have some bad features. I will mention two. 
First when using the Cooper-Frye formula, we sometime 
meet negative terms {da^p^ < 0) corresponding to parti- 
cles re-entering the fluid. However since they presumably 
had stopped interacting (being in the frozen out region) , 
they should not re-enter the fluid and start interacting 
again. 
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FIG. 1: In Cooper-Frye formula, the expression p'^da^ may 
be negative. 

So usually one removes these negative terms from the 
calculations as being unphysical. However by doing this, 
one removes baryon number, energy and momentum from 
the calculation and violates conservation laws. It is not a 
negligible problem, as shown in the figure HJ In the code 
SPHERIO ^ , it can be a 20% overestimate of particle 
number. There are some ways to avoid these violations 
but none is completely satisfying [2ll |. 

The second problem is the following: do particles really 
suddenly stop interacting when they reach a certain hy- 
persurface? Intuitively no, this must happen over a mean 
free path. This is corroborated b y re sults from simula- 
tions of microscopical models [sol Isil Ij^: the shape of 
the region where particles last interacted is generally not 
a sharp surface as assumed for sudden freeze outs. Some 
exceptions might be heavy particles in heavy systems or 
the phase transition hypersurface. 

We postpone the discussion of the second problem to 
section 4 and turn to the first problem. 

2.2 Improved freeze out 

In this section, we adopt the sudden freeze out picture 
and seek ways to incorporate conservation laws [2l|. 



FIG. 2: Rapidity distribution of particles freezing out on the 
27. out = OATo isotherm in the Landau model. The solid line 
corresponds to all the contributions in formula 1, the dashed- 
dotted line represents contributions from space-like parts of 
the isotherm and the dashed line contributions from time-like 
parts |3^. l34l . In this last case, note the negative contribu- 
tions at y=0. 

We suppose that prior to crossing the freeze out sur- 
face (7, particles have a thermalized distribution function 
and we know the baryonic current and energy momentum 
tensor, tiq and Tq"^ . We suppose also that after crossing 
the surface, the distribution function is 

rpo{x,p;dan = fFo{x,p)Qip'' da^) (2) 

The Q function selects among particles which are emitted 
only only those with da^p^ > 0. This equation is solved 
in the rest frame of the gas doing freeze out (RFG) in 
figure 13 . We see that according to the value of v = 
dao/da^j^pQ, a region more or less big in the p± — p^ 
space can be excluded. 
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FIG. 3: The solution to the equation da^p'^ = for v = 
0,65, and -0,35 is given by the left hyperbola, the vertical 
axis and right hyperbola respectivaly. The permitted region 
is localized to the right of the curve in each of these three 
cases. The dashed lines correspond to the case of massless 
particles. 
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We do not know what the shape of fFo{x,p) is. To 
simpHfy, we first suppose that 

fFo{x,p) = j^e^p y — f — J (3) 

with — 7(1, V, 0, 0) and fi is the baryonic potential. 
This does not mean that fpo is thermahzed but simply 
that we choose a parametrization of the thermalized type. 
This parametrization is arbitrary, we discuss later how to 
improve our ansatz. For the moment we use it to illus- 
trate how to proceed in order not to violate conservation 
laws when using the Cooper-Fryc formula. 

It is possible to find expressions for the baryonic cur- 
rent, energy momentum tensor and entropy current cor- 
responding to |2] in terms of Bessel-like functions and for 
massless particules, even analytical expressions as func- 
tion of w, T e /X 

To determine the parameters v,T and /i for matter on 
the post-freeze out side of cr, we need to solve the con- 
servation equations 

[N^da,,] = [T'^^da^] = [T^^da,,] = 0, 
as function of quantities for matter on the pre-freeze out 
side, vq, Tq e fiQ. This being done, we still need to check 
that 

[5^da,]>0ori?=|3^>l 
i.e. entropy can only increase when crossing a. Gener- 
ally, these equations need to be solved numerically but 
for massless particles, they have an analytical solution 
[2l| . For illustration, we show this solution for the case 
of a plasma with an MIT bag equation of state in figure^l 
An interesting result can be seen on the top figure. Nor- 
mally when using a Cooper-Frye formula, the velocity of 
matter pre and post freeze out is assumed to be the same. 
However in the figure, one sees that when imposing con- 
servation laws, matter may be acelerated in a subtantial 
way. For example in case a), vq = 0.2 implies Vfiow = 0.4 
and V = 0.6. In term of effective temperature, there was 
an increase of 60%. 

This example illustrates the importance of taking into 
account conservation laws when crossing a. However, the 
choice of fpo as being parametrized in the same way as 
a thermalized distribution is arbitrary as we mentioned 
already. So we now study a more physical way of com- 
puting this function. 

Consider an infinite tube with the a; < part filled with 
matter and the a; > part empty. At t=0, we remove 
the partition at a; = and matter expands in vacuum. 
Suppose we remove the particles on the right hand side 
and put them back on the left hand side continuously 
so as to get a stationary flow, with a rarefaction wave 
propagating to the left of the matter. 

In the spirit of the continuous emission model pre- 
sented below, the distribution function of matter has two 
components, ffree and Suppose that ffree{x = 

0,p) = and fint{x = 0,p) = ftherm{{x = 0,p). A 



-1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



0.9 




To [MeV] 



FIG. 4: Solution of conservation laws in the case of a plasma. 
Top: V as function of vo (solid line) for a) no = 
1.2/m~^^o = 60AfeV",As = B^^^ = 225 MeV, b) 
no = 0.1/m-^^o = 60 MeV, As = 80 MeV, c) no = 
1.2 fm~^, To = 60 MeV, As = OMeV. (Dashed lines: ve- 
locity of pos- freeze out baryonic flow). 

Middle: baryonic density n as function of no for vq = 0.5, 
To = 50 MeV, a) Ab = 80 MeV (continuous line), b) As = 
120 MeV (dashed-dotted line), c) As = 160 MeV (dashed 
line) . 

Bottom : R, ratio of entropy currents for post and pre- 
freeze matter as function of Tq for a) no = 0.1/m~^,i;o = 
0.5 MeV, As = 80 MeV (solid line), b) no = 0.5/m-^^;o = 
0.5 MeV, As = 80 MeV (dashed line) , c) no = 1.2/m"^,wo = 
0.5 MeV, As = 225 MeV (dotted line). 



4 



simple model for the fluid evolution is 

cos f ,T 



dxfint{x,p)dx 



A 



-fint{x,p)dx, 



da:ffree{x,p)dx = +9 (p^rfg^ ) '^"^ f^nt{x,p)dx. (4) 

where cosOp = ^ in the rest frame of the rarefaction 
wave. A solution for these equations is 

fintix,p) = ftherm{x = 0,p) CXp 



(5) 



and 



1 — exp 



A 



ftherm{x = 0,p) - fint{x,p) 



(6) 



We see that f free tends to the cut thermalized distri- 
bution we saw above when x — > oo. In this model, the 
particle density does not change with x but particles with 
p'^do-^ > pass gradually from fint to ffree- 
To improve this model, we consider 

dxftnt{x,p)dx = -Q{p'^d(Tfj,)^^^^^^^f.M{x,p)dx 

+ [feq{x,p) - fint{x,p)]^dx, (7) 



COS (/ ~* 

dxffree{x,p)dx = +Q{p^da^j) — Y^fint{x,p)dx. (8) 

The additional term in fint includes the tendency for this 
function to tend towards an equilibrium function due to 
collisions with a relaxation distance A'. Due to the loss 
of energy, momentum and particle number, f^q is not the 
initial thermalized function but its parameters neq{x), 
Teq{x) and u^^(a;) can be determined using conservation 
laws. In the case of immediate re-thermalization [fint — 
/eg), for a gas of massless particles with zero net baryon 
number, the solution is shown in figure El One sees that 
the solution //ree is not a thermalized type cut function. 

Even more importantly, this distribution exibits a cur- 
vature which reminds the data on p^ distribution for 
pions. Other explaination for this curvature are trans- 
verse expansion or resonance decays. On the basis of our 
work, it is difficult to trust totally analyses which extracts 
thermal freeze out temperatures and fluid velocities using 
only transverse expansion and resonance decays. 

More details and improvement on how to compute the 
post freeze out distribution function were presented in 
various papers 

Finally, it is interesting to note that models that com- 
bine hydrodynamics with a cascade code also suffer from 
the problem of non-conservation of energy-momentum 
and charge or inconsistency 0| , as discussed in ■ 
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FIG. 5: ffree{Px) ( equivalent to the fpQ in the previous 
section), computed at py — 0,x — lOOA, To = 130 MeV. 



3. SEPARATE FREEZE OUTS 

In this section, we suppose that the sudden freeze out 
picture can be used and see how well it describes data. 

3.1 Chemical freeze out 

Strangeness production plays a special part in ultra- 
relativistic nuclear collisions since its increase might be 
evidence for the creation of quark gluon plasma. Many 
experiments therefore collect information on strangeness 
production. 

We can consider for example the results obtained by 
the CERN collaborations WA85 (colhsion S-HW), WA 
94 (collision S+S) and WA97, later on NA57 (colhsion 
Pb+Pb). One can combine various ratios to obtain a 
window for the freeze out conditions compatible with all 
these data. The basic idea is simple, for example: 

A/A ~ exp (^'5' _ ^^^^ value. (9) 



T \f.out 



(neglecting decays.) 



In principle, this equation depends on three variables. 
However, supposing that strangeness is locally conserved, 
this leads to a relation fj.s{fJ.b-, T), then given a minimum 
and a maximum values, the above equation gives a re- 
lation Tf,out{pB f.out)- We show typical results in fig- 
ure El The parameter 7 in this figure is basically a phe- 
nomcnological factor, which indicates how far from chem- 
ical equilibrium we are, it is introduced in front of the 
factors 6*^=/-^ where /i^ is the s quark chemical poten- 
tial, /is = /is/3 — /i5. (There exists a study by C.Slotta 
et al. |l43| motivating this way of including 7.) It can 
be seen that if 7 = 0, 7 it is possible to reproduce the 
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FIG. 7: Phase diagram with lines indicating values ol Ireeze 
out parameters at various energies obtained from particle 
abundances (solid line) and particle spectra (dashed line) |5^ . 



FIG. 6: Search of a window |44l| of values of Tf.out e pib, /.out 
reproducing WA85 data. This window only exists for 7 < 1 
(for this reference). 



3.2 Thermal freeze out 



experimental ratios A/A, ^/A, ^/A choosing T/.o^t 
and /if,_ /.out in a certain window. This window is located 
around T/.o„t - 180 - 200 MeV, /ib, j.out 200 - 300 
MeV. One can be surprised by such high values since 
particle densities are still high. However these results 
are rather typical as can be checked in the tabic 

TABLE I: Summary of freeze out values from J.SolIfrank |45|. 



collision 


T (MeV) 


lib (MeV) 7. ref. 


S+S 


170 


257 


1 [46] 




197±29 


267±21 


1.00±0.21 [47] 




185 


301 


1 [48] 




192±15 


222±10 


1 [49] 




182±9 


226±13 


0.73±0.04 [50] 




202±13 


259±15 


0.84±0.07 [51] 


S+Ag 


191±17 


279±33 


1 [42] 




180.0±3.2 


238±12 


0.83±0.07 [50] 




185±8 


244±14 


0.82±0.07 [51] 


S+Pb 


172±16 


292±42 


1 [52] 


s+w 


190±10 


240±40 


0.7 [44] 




190 


223±19 


0.68±0.06 [53] 




196±9 


231±18 


1 [54] 


S+Au(W,Pb) 165±5 


175±5 


1 [49] 




160 


171 


1 [55] 




160.2±3.5 


158±4 


0.66±0.04 [51] 



The explaination usually given nowadays is that at 
these temperatures, particles are doing a chemical freeze 
out, they stop having inelastic collisions and their abun- 
dances are frozen. 



Transverse mass distributions obtained experimentally 
(when plotted logarithmically) exhibit large inverse incli- 
nations. These are called effective temperatures. 

In the case of hydrodynamics, these effective temper- 
atures are thought to be due to the convolution of the 
fluid temperature with its transverse velocity, both at 
freeze out. So in particular the effective temperatures 
are higher than the fluid temperature. In addition, the 
effective temperature should be larger the larger the par- 
ticle mass is , since the "kick" received in momentum, 
^ Tnvj''^^'l:, due to transverse expansion, is larger (this 
argument is only valid for the non-rclativistic part of 
the spectrum i.e. p± << m; the effective temperature 
does not depend on mass for the part of the spectrum 
where p± » m (but for that part of the spectra other 
phenomena might be important). In general, given an ex- 
perimental m± spectrum, there exist many pairs of Tf,out 
and W/'"^^ which can reproduce it. 

To remove this ambiguity, we can compare the m±^ 
spectra for various types of particles (e.g. jSg), or for a 
given type of particle, for example pions, combine the fit 
of the spectrum with results on HBT correlations (e.g. 

A compilation for various accelerator energies of 
values for T/.out obtained from particle spectra is pre- 
sented in figure 13 (dashed line). A typical value at SPS 
is Tf out ~ 125 MeV. One can be surprised by the fact 
that these values for Tf^out are lower than those obtained 
from particle ratios. The usual explaination for this is 
that 125 MeV corresponds to a thermal freeze out of par- 
ticles, i.e. when they stop having elastic collisions and 
the shape of their spectra become frozen. This model 
with a chemical freeze out followed by a thermal freeze 
out is called separate freeze out model. Its parameters 
depend on the energy as shown in figure the possible 
decrease of Tf^out with increasing energy is discussed in 
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0, 01 . Some possible problems of this model {ft tem- 
perature and pion abundance are discussed later). For 
the moment let us see how to incorporate this descrip- 
tion in hydrodynamics, since so far it was based only on 
the analysis of two different types of data with thermal 
or semi-analytical hydrodynamical models. 



3.3 Is it quantitatively necessary to modify 
hydrodynamics to incorporate separate freeze outs? 

In [i^l ) we made a preliminary study of whether such a 
separate freeze out model would quantitatively influence 
the hydrodynamical expansion of the fluid. For this, we 
used a simple hydrodynamical model, with longitudinal 
expansion only and longitudinal boost invariance 59] . 

For a single freeze out, the hydrodynamical equations 
are 



de 

dUB 

dt 



P 



t 

T 







The last equation can be solved easily 



nB{to)to 



t 



(10) 



(11) 



Given an equation of state, p{e,nB), we can get e(t) 
and nB{t), solving H1U|I . From them, T{t) and Hsit) 
can be extracted. So, if the freeze out criterion is 
Tf.out — constant, one can see which are the values for 
other quantities at freeze out, for example the values of 
tf.out, MB f.out, ■■■ These values being known, spectra can 
be computed. 

Now let us start again with the previous model, but 
we suppose that when a certain temperature T^h.f.out is 
reached (corresponding to a certain tch.f.out) some abun- 
dances are frozen. To fix ideas, let us suppose that A and 
A are in this situation. In this case, for t > tch.f.out, in 
addition to the hydrodinamic equations above, f llO|l . we 
must introduce separate conservation laws for these two 
types of particles 



dt 



= 

t 

^ = 0. 



dt ' t 

Again it is easy to solve these equations 



nA{t) = 
n^it) = 



n\{to)to 



t 



n]y{to)to 



t 



(12) 
(13) 

(14) 
(15) 



For times t > tch.f.out, we need to solve the hydrody- 
namic equations ( I10|l with an equation of state modified 
to incorporate these conserved abundances. 



We suppose that the fluid is a gas of non-interacting 
resonances. 



27r2 
27r2 



V(T)"+'^ K^inmjT) (16) 



Y.{TT^'^^[m,{nm,lT) (17) 



+ "^K.inm^/T)] 



27r2 



Pi 



^nti.i/T 



Y.{Tr+'^^K,{nmjT) (18) 



n=l 



where nii is the particle mass, gi, its degeneracy and /ii, 
its chemical potential, the minus sign holds for fermions 
and plus for bosons. In principle each particle species i 
making early chemical freeze-out has a chemical potential 
associated to it; this potential controls the conservation 
of the number of particles of type i. For particle species 
not making early chemical freeze-out, the chemical poten- 
tial is of the usual type, /i^ = BifiB+SiUs, where i^J■s) 
ensures the conservation of baryon number (strangeness) 
and Bi (Si) is the baryon (strangeness) number of parti- 
cle of type i. So the modified equation of state depends 
not only on T and fiB but also /iA, ^x,etc. (the nota- 
tion "etc" stands for all the other particles making early 
chemical freeze-out). This complicates the hydrodynam- 
ical problem, however we can note the following. 

li nii — fii » T (the density of type i particle is low) 
and nii >> T, (these relations should hold for all parti- 
cles except pions and we checked them for various times 
and particle types) the following approximations can be 
used 



9t 
27r2 



miT) 



3/2(f,i-mi)/T 



Pi 



y.\l + 

nim^ (l + 



15T 105T2 



8m, 128to2 ^ 
3T 15r2 



2m.j 8m/ 



(19) 

(20) 
(21) 



We note that Ci and pi are written in term of 
and T. Therefore we can work with the variables 
T, ^BT'nf^,nj^,etc, rather than T, /.t^, /iA, /^aj The 
time dependence of riA, Ha; is known as discussed al- 
ready. So the modified equation of state can be computed 
from t, T and ^ b ■ 

In figure |S1 we compare the behavior of T and jiB 
function of t, obtained from the hydrodynamical equa- 
tions using the modified equation of state and the unmod- 
ified one. We see that if the chemical and thermal freeze- 
out temperatures are very different or if many particle 
species make an early chemical freeze out, the thermal 
freeze out time is quite affected. Therefore it is impor- 
tant to take into account the effect of the early chemical 
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FIGURE 2 



t [fm] 

FIG. 8: Hb and T as function of time in the case where all 
particles have simultaneous freeze-outs (dashed line) and (I) 
all strange particles in basic multiplets make an early chemical 
freeze-out (continuous line), (II) all strange particles except 
K and K*'s make an early chemical freeze-out (dotted line). 



freeze-out on the equation of state to make predictions 
for observables which depend on thermal freeze-out vol- 
umes (which are related to the thermal freeze out time), 
for example particle abundances and eventually particle 
correlations. 

If the chemical and thermal freeze-out temperatures 
are not very different or if few particle species make the 
early freeze out, one can proceed as follows. One can 
use an unmodified equation of state in a hydrodynami- 
cal code and to account for early chemical freeze-out of 
species i, when the number of type i particles was fixed, 
use a modified Cooper- Frye formula 



da^p^'f{x,p). (22) 



The second factor on the right hand side is the usual one 
and it gives the shape of the spectrum at thermal freeze- 
out, the first factor is a normalizing term introduced such 
that upon integration on momentum p, the number of 
particles of type i is Nji Tnh. f)- As an illustration, using 
HYLANDER-PLUS [13, we show in figure that both 
the shapes of rn^ spectra and abundances can be repro- 
duced for Tch.f. = 176 MeV and Tth.f. = 139 MeV, while 
simultaneous freeze-outs at Tchj. — Tth.f. — 139 MeV 
would yield the correct shapes but too few particles. 

Hirano and Tsuda |2^ confirmed the importance of 
including separate freeze outs in hydrodynamical codes. 
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FIG. 9: m_L spectra: NA49 data and HYLANDER-PLUS 
results. Dashed-dotted curves: single freeze out at Tf.out ~ 
139MeV. Solid curves (resp. dashed): chemical freeze out at 
Tf.out. ch. = 176 (resp. 184) MeV and thermal freeze out at 
Tf.out. th. = 139 MeV (see text). 
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in particular they studied elliptic flow and HBT radii at To compute u'^, ^ and T, we must solve the equations 

RHIC. This is also consistent with results in (61| . of conservation of energy-momentum and baryonic num- 



4. CONTINUOUS EMISSION 
4.1 Formalism and modified fluid evolution 



be 



r 



D^Tt"" = (28) 
D^(nfcU^) = 0. (29) 



In this section, we present a possible way out of the 
second problem mentioned above. In colaboration with 
Y.Hama and T.Kodama, I made a description of particle 
emission jH^, which incorporates the fact that they, 
at each instant and each location, have a certain proba- 
bility to escape without collision from the dense medium 
(said in the same terms as above: there exists a region in 
spacetime for the last collisions of each particle). So the 
distribution function of the system in expansion has two 
terms //ree, representing particles that made their last 
collision already, and fmt , corresponding to the particles 
that are still interacting 

f{x,p) = ffree{x,p) + fint{x,p) (23) 

The formula for the free particle spectra is given by 

Ed^N/dp^ = J d''xDf,[p''ffreeix,p)] (24) 

(neglecting particles that are initially free; note that if it 
were not the case, the use of hydrodynamics would not be 
possible). Df^ indicates a four-divergence in general co- 
ordinates. This integral must be evaluated for the whole 
spacetime occupied by the fluid. 

This way, we see that the spectra contains information 
about the whole fluid history and not just the time when 
it is very diluted. (This formula reduces to the Cooper- 
Frye formula in an adequate limit). 

We can write 

ffree=Vf = r/{l-V)hnt. (25) 

V{x,p) = ffree/ f, the fraction of free particles, can be 
identified with the probability that a particle of momen- 
tum p escapes from x without collisions, so to compute 
this quantity we use the Glauber formula 

'P{x,p) — exp n{x')(TVreidt'^ ■ (26) 

We suppose also that all interacting particles are ther- 
malized, so 

ftnt{x,p) = fth{x,p) = g/{2'Kf 

X l/{exp[(p.u(x) - iJi{x))lT{x)] ± 1},(27) 

where u*^ is the fluid velocity, /i its chemical potential 
(/i = ^bB + ^sS with B the baryonic number of the 
hadron, 5, its strangeness and /is and /is the associated 
chemical potentials) and T its temperature at x. 



In the figure [TUI examples of solution are given. We 
can see that the fiuid evolution with continuous emis- 
sion is different from the usual case without continuous 
emission. For example, and as expected, the temperature 
decreases faster since free particles carry with them part 
of the energy-momentum. 

In principle we have all the ingredients to compute 
(|24|l . However there exist two problems: 1) numerically, 
in the equation 1251 we can have divergencies if V goes 
faster to 1 than /j„t goes to zero 2) the hypothesis that 
fint is termalized (cf. ea. (|27|) ') must loose its validity 
when V goes to 1. To avoid this problem, we divide space- 
time in eq. 124(1 . in two regions: the first with V > Pj^ 
and the second with V < Vj^, for some reasonable value 
of Pjr. Using Gauss theorem, the second part reduces to 
an integral over the surface V = Vj^ (which depends on 
the particle momentum) 

/l = / da^p^ffree (30) 

= I da^p^ftH (31) 

There is still a certain fraction of interacting particles, 1 — 
T'jF of the total, on this surface, it is these particles which 
in principle turn free in the region V > Vj^. To count 
them, we suppose that they are rather diluted (i.e. Vjr 
is rather large) and we can apply a Cooper-Frye formula 
for them 

h ! da^pf'f.^t (32) 
= [ da^p^fth (33) 

So finally 

E-r^ ^h+h-^ —-^ / da^p^^fth. (34) 

dp"^ l-Vjr Jv=V:r 

It is this formula which is used below, with Vy^ — 0.5 
(but we tested the effect of changing this value) and co- 
ordinates adequate for the geometry of the problem. It is 
similar to a Cooper-Frye formula however one must 
note that the condition V = Vjr depends not only on the 
localization of a particle but its momentum, which as we 
will see has interesting consequences. 
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4.2 Comparison of the continuous emission and 
freeze out scenarios 



In this section, I compare the interpretation of experi- 
mental data in both models. 

a. Strange particle ratios 

We saw above that for the freeze out mechanism, strange 
particle ratios give information about chemical freeze out. 
Now we see how to interpret these data within the contin- 
uous emission sceanrio In this case, 
the only parameters are the initial conditions Tg e 
and a value that we suppose average for 7. We therefore 
fix a set of them, solve the equations of hydrodinamics 
with continuous emission, compute and integrate in 
the spectra given by H34|l for each type of particles (we 
include also the decays of the various types of particles 
in one another). In a way similar to freeze out but now 
with initial values instead of freeze out values, we get in 
figure [TTl a window of initial conditions which permits to 
reproduce the various experimental ratios. (We also look 
at other ratios than those shown in the figure, tested the 
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FIG. 10: Fluid evolution supposing longitudinal boost 
invariencelH^. Top two: for a light projectile such as S, T 
and as function of radius; as a first approximation trans- 
verse expansion is neglected. ( Solid lines correspond to a 
fluid with continuous emission, dashed line lines to a fluid 
without continuous emission, i.e. the usual case.) Bottom 
three: for a heavy projectile such as Pb, T, /is and fluid ve- 
locity as function of radius; transverse expansion is included 
(times: 1,4,7,10 fm). 



effect of changes in the equation of state, cross section, 
initial time, type of experimental cutoff.) 

We see therefore that the initial conditions necessary 
to reproduce the WA85 data are 



(35) 



These values may seem high for the existence of a 
hadronic phase, lattice gauge QCD simulations seem to 
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FIG. 11: Window in initial conditions allowing to reproduce 
WA85 data, for an equation of state with excluded volume 
corrections |66(). 

indicate values smaller for the quark-hadron transition. 
Here we can note that 1) values of QCD on the lattice 
are still evolving (problems exist to incorporate quarks 
with intermediate mass, include /if, ^ 0, etc. 2) Our own 
model is still being improved and we know for example 
that the equation of state affects the localization of the 
window in initial conditions compatible with data. 

In figure El the same kind of analysis was done for 
the data WA97. With the reservation in the caption of 
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FIG. 12: Window in initial conditions reproducing WA97 
data in the case of an equation of state without excluded 
volume corrections. With volume corrections, the window is 
lower (gQj. 

the figure, we see that the initial conditions are not very 
different from the one above. 



We can therefore conclude that the interpretation of 
data on particle ratios lead to totally different infor- 
mation according to the emission model used. For the 
freeze out model, we get information about chemical 
freeze out while for continuous emission, we learn about 
the initial conditions. 

b. Transverse mass spectra 

For freeze out, we saw that these spectra tell us about 
thermal freeze out (temperature and fluid velocity). Now 
for continuous emission let us see how to interpret these 
same spectra In this case, the initial 

conditions were already determined from strange parti- 
cle ratios, they cannot be changed and must be used to 
compute the spectra as well. For example, in figure [T^ 
various spectra are shown, assuming Tq — ^ibo ~ 200 
MeV, and compared with experimental data. 

This comparison should not be considered as a fit but 
as a test of the possibility of interpreting various types of 
data in a self-consistent way with continuous emission (in 
particular note that our calculations were done assuming 
longitudinal boost invariance). 

We learn various information from this comparison. 
In these figures, we do not take into account transverse 
expansion. With the S-l-S NA35 data, we note that 
the heavy particles and high transverse momentum pi- 
ous have similar inverse inclinations Tq ~ 200 MeV. The 
particles heavier than pions, due to termal suppression, 
are mainly emitted early when the temperature is ~ To- 
For lower temperatures, there are still emitted (and more 
easily due to matter dilution) but their densities are quite 
smaller (this is what is called thermal suppression) and 
their contributions as well. The high transverse momen- 
tum pions have large velocity and (if not too far away 
from the outer fluid surface) escape without collision ear- 
lier than pions at the same place but with smaller veloc- 
ity, so these high transverse momentum pions also escape 
at ~ Tg. Pions are small mass particles and are little af- 
fected by thermal suppression. This way, they can escape 
in significative number at various temperatures. This is 
reflected by their spectra, precisely its curvature. (In our 
calculation, decays into pions are not included, this would 
flU the small transverse momentum region and improve 
the agreement with experimental data.) 

The S+S WA94 data also indicate that continuous 
emission is compatible with data. Finally, the S-l-W 
WA85 data seem to indicate that perhaps somewhat dif- 
ferent initial conditions or a little of transverse expansion 
might be necessary to reproduce data. 

Our calculations including transverse expansion 
indicate that little transverse expansion is compatible 
with data for light projectile. This is understandable 
noting that the effective temperature of spectra are 
already of order the initial temperature ~ 200 MeV. 
In the case of heavy projectile, the situation is dif- 
ferent. The various types of particles have different 
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temperatures, all well above 200 MeV. In this case, 
we must include transverse expansion to get consis- 
tency with data. An example is shown in figure H"!] with 
the same values of Tq and ^ibfi than previously, 200 MeV. 
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FIG. 13: Using longitudinal boost invariance and no trans- 
verse expansion, we compare our predictions, top to bottom, 
with NA35 data {S+S, all rapidity), with WA94 data (S-fS, 
midrapidity, large transverse momenta) and with WA85 data 
(S-l-W, midrapidity, large transverse momenta) [6^. 



FIG. 14: WA97 data and example of comparison with contin- 
uous emission with transverse expansion (not a least square 
fit) for To = fib.o = 200 MeV 



c. The case of the 

The effective temperatures in the case of heavy projec- 
tiles, particularly for the 51, have attracted a lot of atten- 
tion for the following reason. In the usual hydrodynamics 
with freeze out, as we saw, it is expected that effective 
temperatures increase with mass. In this context it was 
difficult to understand why the f2 had an effective tem- 
perature much lower than other particles in figure 1151 
A possible explaination within hydrodynamics with sep- 
arate freeze outs, is that the made its chemical and 
thermal freeze out together, early. Van Hecke et al. 
argued that this is reasonable since it is expected that 
il has a small cross section (because there is no chan- 
nel for fliT resonance — > riir) and showed that the 
microscopical model RQMD can reproduce the data. 

In our model originally we had used the same value of 
the cross sections to compute the escape probability for 
the various types of particles. This is not expected and 
indeed in this case, continuous emission, does not lead 
to good predictions as shown in figure [T7k . Therefore in 
the spirit of microscopic models, we also show our pre- 
dictions in figure 117b . for continuous emission and the 
following cross sections: ct^tt crvrei >Tnv'^ 1/"^-^, 

cttytt crvrei >N7r^ 3/2 < aVrei >7r7r (usiug addi- 
tive quark model estimate), ctatt ~< (JVrei >A7r~ lj2 < 
aVrei >7r7r (usiug additive quark model estimate [70l |) and 
apj ~< (TVrei >n7T^ 1/2 < cTVrei >n-k (usiug estimate in 
|7l|). The predictions now are in agreement with data. 
However, the cross sections are poorly known and our 
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FIG. 15: Compilation of experimental data on effective tem- 
peratures in the case of heavy projectile and predictions for 
usual hydrodynamics |6^ . (The various experiments have dif- 
ferent rapidity and transverse momentum cutoffs). In general, 
deuteron is left outside the hydrodynamical analysis: it is lit- 
tle bounded and should form later in the fluid evolution from 
coalescence of a neutron and a proton with similar moment. 
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FIG. 16: Predictions from RQMD [69l | in the same experi- 
mental situation as the previous figure. 



results are sensitive to their values. 

Recently, with new data by NA49, WA97 ^ and 
NA49 [2J| came to the conclusion that in fact there 
is no need for early joint chemical and thermal freeze 
outs for the fi: all their spectra can be fitted with a 
simple hydro inspired model as seen in figure 1181 The 
previous difficulty for WA97 came from the fact [t^ 
that SI was observed at high p±. However now, STAR 
has problems 76] to fit with a simple hydro inspired 
model the 5 together with 7r,p, if, A and would need to 
assume early joint chemical and thermal freeze outs for 
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FIG. 17: Comparison of experimental data on effective tem- 
peratures in the case of a heavy projectile with predictions 
from continuous emission. Top: all cross section supposed 
equal. Bottom: more realistic cross section values (cf. text) 

mill. 



the 2, as shown in figure El More recently, it has been 
noted 01 by NA49, that their conclusion depends on 
the parametrization used and in |7^ , NA57 argues that 
due to low statistics, it is not clear what conclusion can 
be drawn for the fl. In |79| . an attempt was made to 
reproduce with a single thermal freeze out temperature 
in a hydrodynamical code, all transverse mass spectra 
at a given energy. 

d. Pion abundances 

We showed that strange particle ratios can be reproduced 
by a model with chemical freeze out around 180 MeV. 
The abundances too can be reproduced. The problem is 
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FIG. 18: Data and hydro inspired fit for all particles including 
top, NA49 and bottom, NA57 US, III- 



that the pion number is too low. This was noticed by 
Davison et al. [s^l, as shown in their table HTl reproduced 
below 



TABLE II: Comparison between NA35 data and previsions 
from a thermal model 1801. 



Th. 

Exp. 



K°s K+ K- A 



A 



10.7 14.2 7.15 8.2 1.5 23.2 56.9 
10.7 12.5 6.9 8.2 1.5 22.0 92.7 
±2.0 ±0.4 ±0.4 ±0.9 ±0.4 ±2.5 ±4.5 



There are various ways to try to solve this problem. 

1. It can be argued (in a spirit similar to Cleymans 
et al. [U) that strange particles do their chemical 
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FIG. 19: A hydro fit to transverse mass spectra leads to high 
Tf.out for H and lower for tt,p,K,A [t^ . 



freeze our early around 180 MeV and their thermal 
freeze out around 140 MeV but the pions do their 
chemical and thermal freeze outs together around 
140 MeV. In fact Ster et al. js^ manage to repro- 
duce data by NA49, NA44 and WA98 on spectra 
{including normalization i.e. abundances) for nega- 
tives, pions, kaons and protons and HBT radii (see 
next section) with temperature around 140 MeV 
and null chemical potential for the pions. On the 
other side, Tomasik et al. [s^ say that for the same 
kind of objective, they need temperatures of order 
100 MeV and non zero chemical potential for the 
pions. So it is not clear if to reproduce the pion 
abundance, it is necessary to modify the hydrody- 
namics with separate freeze out, including pions out 
of chemical equilibrium or not. 

2. Gorenstein and colaborators [s^ Issl Is^ studied 
modifications of the equation of state, precisely 
they included a smaller radius for the volume cor- 
rections of the pion. 



3. Letessier et al. made a serie of papers |53, |SV| ar- 
guing that the large pion abundance is indicative 
of the formation of a quark gluon plasma hadroniz- 
ing suddenly, with both strange and non-strange 
quarks out of chemical equilibrium js^ ■ 

Given the difficulty that freeze out models have with 
pions, it is interesting to compute abundances with con- 
tinuous emission models. In the table, results and NA35 
data from S±S are shown (data are selected at midrapid- 

ity). 

It can be seen that in the continuous emission model, 
a larger number of pions is created. In (SQj . we related 
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TABLE III: Comparison between experimental data, results 
for continuous emission at Tq — /^(,.o — 200 MeV and freeze 
out at To — fj,b.o ~ Tf.out = fib f. out — 200 MeV. 



and 





experimental value 


continuous emission 


freeze out 


A 


1.26±0.22 


0.96 


0.92 


A 


0.44±0.16 


0.29 


0.46 


p-p 


3.2±1.0 


3.12 


1.32 


h- 


26±1 


27 


15.7 


K% 


1.3±0.22 


1.23 


1.06 



this increase to the fact that entropy increases during the 
fluid expansion. This is due to the continuous process of 
separation into free and interacting components as weU 
as continuous re-termaUzation of the fluid. In contrast, 
in the usual hydrodynamic model, entropy is conserved 
and related to the number of pions; so in this usual 
model, observation of a large number of pions implies a 
large initial entropy and is indicative of a quark gluon 
plasma lUE^, 88]. In our case, a large number of pions 
does not imply a large initial entropy and the existence 
of a plasma. It can be noted that the interpretation 
of data is quite influenced by the choice of the particle 
emission model. 

e. HBT 

Interferometry is a tool which permits extracting infor- 
mation on the spacetime structure of the particle emis- 
sion source and is sensitive to the underlying dynamics. 
Since pion emission is different in freeze out and contin- 
uous emission models, it is interesting to cornpare their 
interferometry predictions (for a review see [90j). This 
was done in ref. [9l|. 

In this work, the formalism of continuous emission 
[b^ IH^j l was extended to the computation of correlation 
functions. Precisely, we computed 



C{ki,k2) = C{q,K) = l + 



\G{q,K)\ 



G(fcl,fcl)G(fc2,fc2) 



(36) 



where qt" = k'^; - e Rt' = i(A:f -I- fc^). 

In the case of freeze out (in the Bjorken model with 
pseudo-temperature[9^ ) 

G(fci,fe) = 2 < ^ > {^^J,{qTRj.)}Ko{0 (37) 



dy qrRr 



G(fc„fc.) = 2^Xo(^). 



(39) 



In the case of continuous emission in the Bjorken model 
with pseudo-temperature 



G{q,K) = 



1 



(2^)3(1- 7^^) 7o 

•Rt 



2-K p + oo 

d(f) / drj 



{ p dp Tjr Mt cosh(y — rf) 
Jo 



X e 



i['<'3^{qo cosh77-g_L sinh)))— pgr 003(0—0,)] 



X e 



X e 



TdrpyrKx cos (f) 

i[T{qo coshr(-gi, sinh r;) — pjfiJt cos(0— 
-Mr coBh{Y-r]) /Tp^{x) 



(40) 



with Mt = y/K'^ + AP, Kt = \{ki + fe2)T , = 
Kfj^Rf^ — — jqfj,q^, Y is the rapidity corresponding 
to K, (j) is the azimuthal angle in relation to the direc- 
tion of K and (pq is the angle between the directions of 
q and K. r-p^ and p-p^ are determined hy V = Vj^. 
Tps{x) = l,42r(a;) - 12, 7 MeV. 

In [91|, a few idealized cases were studied and then 
some cases more representative of the experimental situ- 
ation, were presented. For example, instead of C{q,K), 
we computed 

{C{qL)) = 1 + 

rZo c '^^^ c c <^M,)\GiK,,)\^ 



/_™^ d/fi yj dKT dqs y" dq,C{K,q)G{k,M)Gik2M) ' 



(41) 



(This corresponds to the experimental cuts of NA35). qp, 
qs ^ qL are defined in figure EUl 









'''^J>^q = (qo'qs'q,) - 








z 

1 = "lona" 



where 



'2T 



{miT + m2T) - iT{miT - m2T)]^ + 



2{^+T^)m,Tm2T[co8h{Ay)-l] , (38) 
= 2/1 ~ 2/2, < > indicates average over particles 1 and 

2 



FIG. 20: By definition : Oz is chosen along the beam and 
K in the x ~ z plane. The L (longitudinal) component of a 
vector is its z component, the O ("outward") its x component 
and S ("sidewards") its y component [o^. 
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In a first comparison, we used similar initial conditions 
than above, Tq = 200 MeV (S+S collisions) for both con- 
tinuous emission and freeze out. Results are presented as 
function of QljQo ^tnd 95 in figure f?n In a second com- 
parison, given a curve obtained for continuous emission, 
we try to find a similar curve obtained with the standard 
value Tf,out = 140 MeV, varying the initial temperature 
T(( °"*. The results are shown as function of qL.qo e qs 
in figure 122 From these two sets of figures, it can be 
seen that there are many differences in the correlations 
between both models: shape, heigth, etc. If the initial 
conditions are the same, the correlations are very differ- 
ent. Even more interesting if trying to approximate with 
a freeze out at Tf,out = 140 MeV, the continuous emis- 
sion correlations, it is necessary to assume a Tq very high, 
where the notion of hadronic gas looses its validity. So, 
we can see that if continuous emission is the correct de- 
scription for experimental data, it will be more difficult 
to attain the quark gluon plasma than it looks using the 
freeze out model. So again we conclude that what we 
learn about the hot dense matter created depends on the 
emission model used. ( Note that to actually compare 
with data, transverse expansion has to be included. ) 

More recently, we addressed the "HBT puzzle" at 
RHIC using NeXSPheRIO ^ (for a review on NeX- 
SPheRIO, see jo^). We showed that the use of varying 
initial conditions leads to smaller radii. In addition, 
though continuous emission is not easy to introduce 
in hydrodynamical codes (because H26|l depends on 
the future evolution of the fluid), it was introduced 
in an approximate way and shown to be important to 
reproduce the momentum dependence of the radii. 

f. Plasma 

In all the previous analysis, we assumed that the fiuid was 
initially composed of hadrons with some initial conditions 
Tq, /Xb_o- Given the possibility that a quark gluon plasma 
might have been created already, we must discuss the 
extension of our model to the case were a plasma might 
have been formed. 

Initially, it might be expected that continuous emission 
by a plasma is impeded for two reasons: 1) a hadron 
emitted by the plasma surface, in the outward direction, 
must cross all the hadronic matter around the plasma 
core , so probably it will suffer collisions and when finally 
emitted, it will be emitted by the hadronic gas in the 
way we have considered so far 2) the plasma is formed 
by colored objects that must recombine in a color singlet 
at the plasma surface to be emitted, this makes plasma 
emission more difficult than hadron gas emission. 

For simplicity let us consider a second order phase 
transition. We use for the hadron gas a resonance gas 
equation of state. For the plasma, we use a MIT bag 
equation of state where the value of the bag constant and 
transition temperature are adjusted to get a second or- 
der phase transtion. We get ~ 220 MeV and B ~ 580 
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FIG. 21: Comparison of continuous emission with freeze out, 
when both have the same initial conditions lOlll. 
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FIG. 22: Given a curve obtained for continuous emission, we 
look for initial conditions for freeze out at Tf,out ~ 140 MeV 
leading to a similar curve (oH- 



MeV /m-3. 

Then we solve the hydrodynamics equation (without 
continuous emission as a first approximation) to know 
the localization and evolution of the plasma core. This is 
shown in figure ESI These equations were also solved for 
the case of a hadronic gas for comparison in figure EH 




10 15 20 
r(fm) 

FIG. 23; Evolution of temperature in the case of a hadron 
gas with quark core and initial boxlike matter distribution, 
top and softer, bottom. 

It can be checked that when there exists a plasma core, 
it is quite close to the outside region. Contrarely to the 
reservation 1) above, hadrons emitted by the plasma sur- 
face might be quite close enough to the outside to escape 
without collisions. 

Now for reservation number 2}, we n ote t hat there exist 
various mechanisms [o^, IH, |9^ Il00| proposed for 
hadron emission by a plasma. To start we can assume as 
Vishcr et al. [ooj that the plasma emits in equilibrium 
with the hadron gas. In this case, the emission formula 
by the plasma core+hadron gas would be 



E- 



d(j)dl][ pdp{p^ f freer). 



/ '^dT{pPffreep)lii^^^(^r)] 
JTn 



(42) 
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FIG. 24; Evolution of temperature in the case of a hadron 
gas and initial boxlike matter distribution, top and softer, 
bottom. 



where Rout(T) is the radius up to which there is matter. 
This is similar to the hadron gas case treated above. The 
difference is in the calculation of V, (which appears in 
ffree), siucc SL liadron entering the plasma core will be 
supposed detroyed. 

In the same spirit as above, a cutoff at 7^ = Pj^ can 
be introduced. Due to the similarity for the spectra for- 
mula with and without plasma core, we do not expect 
very drastic differences if the transition is second order. 
Of course, the case of first order transition must be con- 
sidered (though results from lattice QCD on the lattice 
do not favour strong first order transition). (Note that 
Tc ^ 220 MeV is higher than expected, this might be 
improved e.g. using a better equation of state). 



CONCLUSION 

In this paper, we discussed particle emission in hydro- 
dynamics. Sudden freeze out is the mechanism commonly 
used. We described some of its caveats and ways out. 

First the problem of negative contributions in the 



Cooper-Frye formula was presented. When these con- 
tributions are neglected, they lead to violations of con- 
servation laws. We showed how to avoid this, the main 
difficulty remains to compute the distribution function of 
matter that crossed the freeze out surface. Even models 
combining hydrodynamics with a cascade code have this 
type of problem or related ones [i^. 

Assuming that sudden freeze out does hold true, data 
call for two separate freeze outs. We argue that in this 
case, this must be included in the hydrodynamical code 
as it will infl uence the fluid evolution and the observables. 
Some works |lOlj using paramatrization of the hydrody- 
namical solution suggest that a single freeze out might 
be enough. No hydrodynamical code with simultaneous 
chemical and thermal freeze outs achieves this so far (see 
e.g. our figures EJl. On the other side, (single) explosive 
freez e out is being incorporated in a hydrodynamical code 

Finally, we argued that microscopical models indeed do 
not indicate a sudden freeze out but a continuous emis- 
sion l^^j 1^ 113 ■ We showed how to formulate particle 
emission in hydrodynamics for this case and discussed 
extensively its consistency with data. We pointed out 
in various cases that the interpretation of data is quite 
infiuenced by the choice of the particle emission mecha- 
nisms. The formalism that we presented for continuous 
emission needs improvements, for example the ansatz of 
immediate rethermalization is not realistic. An example 
of such an attempt is 0| . Finally, it is also necessary to 
think of ways to include continuous emission in hydro- 
dynamics. This is not trivial because the probability to 
escape depends on the future. In 94], such an idea was 
applied to the "HBT puzzle" at RHIC with promising 
results. 
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